Giant radiation heat transfer through the micron gaps 
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Near-field heat transfer between two closely spaced radiating media can exceed in orders radiation 
through the interface of a single black body. This effect is caused by exponentially decaying (evanes- 
cent) waves which form the photon tunnel between two transparent boundaries. However, in the 
mid-infrared range it holds when the gap between two media is as small as few tens of nanometers. 
We propose a new paradigm of the radiation heat transfer which makes possible the strong photon 
tunneling for micron thick gaps. For it the air gap between two media should be modified, so that 
evanescent waves are transformed inside it into propagating ones. This modification is achievable 
using a metamaterial so that the direct thermal conductance through the metamaterial is practically 
absent and the photovoltaic conversion of the transferred heat is not altered by the metamaterial. 

PACS numbers: 78.67.Ch,77.84.Lf,41.20. Jb 



(N 



C/3 

o 

• 1—1 

> . 

Oh. 

q 
o 



> 
o 
p 

rn 
O 



X 



I. INTRODUCTION 

Over the last decade near-field thermo-photovoltaic 
(NF TPV) systems (e.gii"— ) are often considered in the 
modern literature as a promising tool for the field recu- 
peration from high-temperature sources (industrial waste 
heat, car engine and exhausting pipe, etc.). Also, NF 
TPV systems are used as precise temperature profile sen- 
sors (e.giii) . The operation principle of NF TPV systems 
is based on the use of the evanescent spatial spectrum, 
i.e. on the use of the energy of infrared fields stored at 
nanometer distances from the hot surface. In presence 
of another body in the near vicinity of the hot surface 
the well-known photon tunneling phenomenon arises be- 
tween two surfaces which leads to the dramatic increase 
of the heat transfer compared to the value restricted by 
the back-body limit. This transfer can be enhanced if 
the hot medium, denoted as medium 1 in Fig. [T] (a), 
and the photovoltaic material, denoted in Fig. [T] (a) as 
medium 3, both possess negative permittivityii^i^ii^^. The 
enhancement holds at a frequency where Re(ei_3) ~ —1 
(when medium 2 is free space) and is related to the ex- 
citation of coupled surface-plasmon polaritons (SPP) at 
the interfaces of media 1 and 3^. 

In spite of their relatively high efficiency, NF TPV sys- 
tems suffer strong technological drawbacks which restrict 
their applicability and industrial adaptation. First, it 
is difficult to create flat and clean surfaces of Media 1 
and 3 separated with a nanogap (the roughness should 
be much smaller than the gap thickness d). Second, the 
cryogenic cooling or/and vacuum pumping is required for 
their operation since the thermal phonons suppress the 
photovoltaic conversionl''' and the thermal conduction in 
the air is very high if d is smaller than the mean free path 
of the air molecules (30 — 50 nm) or comparable with it. 

Recently, works on the theory and design of so-called 
micro-gap TPV (MTPV) systems have appeared (see e.g. 
ii>iiiS). MTPV systems in which the gap d between the 
photovoltaic and hot surfaces is smaller than the opera- 
tional wavelengthes but comparable with them occupy an 
intermediate place between conventional (far-field) TPV 
and NF TPV systems. This concerns also their efficiency 
which is much less than that of NF TPV systems though 



much higher than that of far-field TPV systems. The 
lack of efficiency is partially reimbursed by strong tech- 
nological advantages. Attempts to increase the efficiency 
of MTPV systems are knownii^iiS related with the use 
of nanoantennas arrays created on both hot and photo- 
voltaic surfaces or insertion of a photonic crystal layer 
inside the gap, but the enhancement is not as signifi- 
cant as to justify these complications. In principle, one 
can strongly increase the radiation heat transfer through 
the gap partially filling it with a lossless negative-index 
metamateriaiii, however high losses are inherent to such 
metamaterials and for substantial gaps and known meta- 
materials this mechanism is not very efficient. 

In the present paper we suggest a new paradigm for 
MTPV systems which should make such structures more 
competitive than NF TPV systems. We suggest to fill 
in the gap between Media 1 and 3 with the metama- 
terial 2, performing the transformation of the evanes- 
cent spatial spectrum into propagating one. Metamate- 
rials performing this manipulation with electromagnetic 
waves are, e.g., the so-called indefinite medic^ - uniaxial 
materials in which the axial and tangential components 
of the permittivity tensor has different signs. For the 
infrared range such metamaterials can be performed as 
arrays of aligned carbon nanotubes (CNT)i^ii^. Arrays 
of aligned metallic CNTs transform incident p-polarized 
infrared waves (propagating and evanescent) into quasi- 
TEM waves propagating along CNTs with quite small de- 
cay. Unlike the SPP enhancement this effect holds over 
the whole mid IR rang o^^'^'^ . 
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FIG. 1: (Color online) (a) - Illustration to the general problem 
formulation, (b) - Interdigital arrangement of CNT allows to 
get rid of the thermal conductance. 
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Arrays of single-wall aligned carbon nanotubes with 
metallic properties are fabricated by many research 
groups and used as field emitters^-, biosensorai^ and 
antennasi^iii. Such parallel CNTs can be few microme- 
ters long and by two-third of their length free standing 
(one third is submerged into a dielectric matrix)^ ^. It is 
possible to perform medium 2 so that the direct thermal 
conduction through CNTs is avoided. A possible scheme 
is shown in Fig. [1] The area filled with medium 2 is di- 
vided by three layers: d = hi + h2 + h^. Since there is no 
mechanical contact of any CNT with both media 1 and 3 
and, since the array density is assumed to be very small, 
the thermal conductance through the array of CNTs will 
be smaller than the thermal conductance through the air 
gap (for d > 1 fim the last one does not lead to the critical 
suppression of the photovoltaic conversion). 



II. THEORY 

The radiation heat transfer to medium 3 is calculated 
through the ensemble- averaged Poynting vector (S^'^) cre- 
ated by medium 1 at the input of medium 3 (plane z — d) 
minus the backward Poynting vector (S^^) (calculated 
at z = 0)^. The spectral density of the total heat net 
transferred between medium 1 with temperature Ti and 
medium 2 with temperature T2 equals^: 

[{Sf{K,uj,Ti)) - (Sf (A:„w,r2))] k^dk^. 

(1) 

In ([1]) we neglect the heat radiation and absorption in 
medium 2. This approximation is justified for CNT ar- 
rays of small density (ratio diameter/period). 

To calculate the Poynting vector (Sj^) (the value (S^^) 
is obtained similarly) we first find the incident Poynting 
vector {S\{kx, oj)) calculated for the semi-infinite medium 
1 neglecting the reflection from the interface between me- 
dia 1 and 2. Then we apply the transmission matrix 
method which allows us to express (Sj^) through the 
incident Poynting vector (Sj(fcx,tij)) and known param- 
eters of media 2 and 3. Following to^ we start from 
nonhomogeneous Maxwell's equations with random elec- 
tric current density sources j and apply the fluctuation- 
dissipation theorem (see e.g. in^-^) for the ensemble- 
averaged bulk current density: 



(j„(r,w)j,*(r',a.')) = 

^ujeae"{uj)5mn5{^ - v')5{uj - T). 



(2) 



In Eq. ^ e" = Im(ei), jm and j„ {m, n — 1, 2, or 3) are 
X, y ov z component of j, 5,„„ is the Kronecker symbol, 
5{x) is the Dirac delta function, and 0(w, T) is the mean 
energy of the Planck oscillator 

Q("'^)^ exp(;.Jfc'lT)-l 

where h is the reduced Planck constant, fc^ is the Boltz- 
mann constant, and T is the absolute temperature of the 
medium. According to the ergodic hypothesis, the spec- 
tral energy flux of p-polarized waves is expressed as^": 

1 

(Si(r,w))-y^ -{KeE.,{v,uj)H;{v,J))d^' (4) 



To find the ensemble-averaged product {ExH*) and in- 
tegrate it over uj' is easy taking into account the homo- 
geneity of the structure in the {x — y) plane. It allows us 
to present an elementary bulk current source in a form of 
a harmonic current sheet: 



j(z) =jo{z')S{z - z')eJM-fc=^^-'=y?'). 



(5) 



Without loss of generality we can put here ky = 0. Solv- 
ing Maxwell's equations with the source defined by ([U 
we obtain horizontal field components created at {x, z) 
by the source located at z': 



'^^x = '^^Uoxkiz - jozkx)e- 



= w-{-ioxkiz + iQzkx)& 



j[ujt — kxX — k-iz(z — z')] 
j[uit-k^x-ki;,(z-z')] 



where kiz = \J fc^e^ — fc^, r\ = y/eaJJ^- The dependence 
on X disappears in the cross product Q . Taking into ac- 
count that in ^ there is no correlation between different 
current sheets we can present the integrand in ^ as: 

{ExH;) = ^LZ^/°^e.('=--feL)^'e.(--')*([j„^(^,^') 
kiz - iQz{^,z')kx\[-jax{^' ,z')ki,_ + jn^{uj' , z')kx])dz' . 

Further, we can rewrite Eq. ^ in the form comprising 
the current density amplitudes: 



{jo^{z',uj)j*Jz',uj')) = 
|weoe"(w)(5™„(5(w - uj')e{uj, T). 



(8) 



Substituting ([7]) into (|4]) we easily perform integrating 
over Lj' and ensemble averaging using ([SJ. Integrating 
over z' we take into account that 

/ dz'exp [j(fci2 - kl^)z'] = -- (Imfci^) ^ . (9) 

J -oc ^ 

The incident Poynting vector created at the point z — Q 
at the frequency lo by sources of thermal radiation dis- 
tributed in medium 1 yields as follows: 



e'/(a.)/27r .(^2 ^ fc^^fc*j0(^^ j.) + ^^^^ 



(10) 



eifci^Im(fciz) 



Now we have to express the Poynting vector trans- 
mitted from medium 1 to medium 3 taking into ac- 
count the wave reflections from two interfaces. Assuming 
that Maxwell's boundary conditions are satisfied at both 
boundaries we can present it through the wave trans- 
mission coefficient r and transverse wave impedances of 
media: 



(Sf(A:.,..)) = i(Si(fc.,a.))|rpg+c.c., (11) 



where 



-E,x/H,y^7^h,/k (i== 1,2,3). 



(12) 



The coefficient r can be obtained through transmission 
matrix components of medium 2 in a form'^^: 



(13) 



Mil + M12/Z3 + M21Z1 + M22^l/^3 ' 
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where Mmn can be found through the z-component 
of the wave vector in medium 2: 



11 



M22 = cos k2zd, A/12 ~ 3Z2 sin fc2zd, 



M21 = sinfc2^d. 



(14) 



If the medium 2 is free space fc2z = — j y^k^ — and for 
evanescent waves |r| decays exponentially when k^d > 
1 (beyond the region of SPP^'^^). However, medium 2 
formed by CNTs even for very large values oik^/k possess 
real-valued k^ (this is so if we neglect losses in CNTs) 
and the contribution of spatial frequencies corresponding 
(in free space) to the evanescent spectrum is important. 
In our calculations we took into account real losses in 
CNTs, however these losses practically do not alter the 
transmission of waves with large k^- 

The applicability of the approach based on the transfer 
matrix for finite-thickness slabs of aligned CNTs (prac- 
tically the validity of Maxwell's boundary conditions at 
their boundaries) was checked in^. Formula ((T2|) and 
expression (IT5|) for the transfer matrix allow high accu- 
racy when calculating the transmission through a layer 
filled with a low-density array of aligned CNTs. Since the 
backward flux {S'l^{kx,uj)) calculates similarly to (Sj^), 
all we need to complete our model is to relate k2z to 
kx for given lo. For this we used the model of the indefi- 
nite medium with permittivity components developed for 
aligned CNTs in^^. Taking into account the interdigital 
arrangement of CNTs depicted in Fig. [T] the transfer 
matrix B of the whole gap is calculated as the product 
B = Ml X M2 X Ml, where Mi, M2 are transfer matrices 
of layers hi, /12 and /13, respectively. Matrix M2 corre- 
sponds to a twice more dense array than matrices Mi ^3 
do. 

Now let us assume that media 1 and 3 have the 
same permittivity given by the Drude formula for heavily 
doped silicon^^: 



(15) 



where £00 ~ 11-6 is the high-frequency limit value of 
the permittivity^^, and 7 is the scattering rate. The 
plasma frequency and scattering rate are expressed as 
ujp = yJWe^~jl(mFeQ) and 7 = e/ (m*/i, respectively, where 
e is the electron charge, N is the carrier concentration, 
m* is the carrier effective mass, and fi is the mobility. For 
n-type heavily doped Si (namely this material is used in 
NF TPV systems enhanced by coupled SPP— i^) the mo- 
bility expression is given as^^ 



M = Ml + 



Ml 



M2 



1 + {N./CrT 



1 + {Cs/Ne) 

1414cm^/Vs, ^2 



(16) 



Cs = 3.42 X 10^" cm 



56.1 

3 



Here /ii = 68.5 cm^/Vs, ^max = 
cmVVs, a = 9.2 X 10^^ cm-3 
a = 0711, P = 1.98, and N^. is the electron concentration. 

Fig. [2] presents the comparison of transmittances jrp 
calculated in presence of the CNT array and in absence 
(vacuum gap). Calculations were done at the frequency 
corresponding to the wavelength A — 7.5 //m. Carri- 
ers concentration for n-doped silicon was taken — 
10^° cm^'^. The width of the gap between media 1 and 3 
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FIG. 2: (Color online) Transmission coefficients versus k^/k 
for evanescent waves in ttie vacuum gap (blue) and for waves, 
propagating via the CNT array (red). 



was equal in this example d — 500 nm and thicknesses of 
sections were as follows: hi = h^ = 50 nm, ft,2 = 400 nm. 
Periods of CNT arrays 01^3 = 20 nm for sections 1 and 3 
and 02 — 10 nm for the central section 2. Radius of the 
CNT was taken equal r — 0.822 nm and the relaxation 
time Tr = 10-^3 s2^. 

Transmission peaks, exceeding unity, correspond to the 
excitation of coupled SPP at the interfaces of media 1 
and 3 (compare with results oS"— i^^). The most impor- 
tant peak corresponds to the case of the vacuum gap and 
is located at spatial frequencies k^/k = (1.5... 1.7). In 
presence of CNTs this peak is suppressed. However, in- 
stead, multiple peaks of |t| appear in presence of CNTs 
and (what is most important) |t|, which in absence of 
CNTs practically vanishes at k^ > lOfc, in their pres- 
ence does not. High spatial frequencies k^/k — 5 ... 25 
correspond to |t| > 0.1. The limit value k^/k = 50 
corresponds to the condition kxai ~ 7r/3. For spatial 
frequencies exceeding this limit value the conversion of 
evanescent waves into propagating ones disappears (see 
also in^^), and corresponding spatial harmonics practi- 
cally do not transport energy. Therefore calculating the 
transferred heat with formula ([IJ we implement integra- 
tion over the Brillouin zone of the CNT lattice with the 
constant ai = 20 nm. Obviously, the integral over k^ in 
formula ([T]) would diverge if we assume that fc^^ — > 00 still 
correspond to waves propagating along the CNT axes. 

Fig. |3] illustrates the dependence of the spectral density 
of the transferred heat g"(A) (the evident analogue of g") 
for the case, when d = 100 nm, hi = h^ — 10 nm and for 
d = 1000 nm, hi ~ hs — 100 nm. The temperatures of 
the hot and photovoltaic surfaces were in both examples 
picked up as follows: Ti = 250° and T2 = 0. Then the 
maximal heat transfer corresponds to At ~ 9.5 fim and 
(S^i) = 0. For the gap d = 100 nm the insertion of CNTs 
gives a modest overall enhancement about 20%. This is 
so because for the gap d — 100 nm the SPP enhance- 
ment (suppressed by CNTs) is still strong. However, for 
the gap d = 1000 nm, where the SPP enhancement is 
negligible, the total thermal transfer offered by CNTs ex- 
ceeds by more than two orders the heat transfer through 
the vacuum gap. This effect opens a new door for the 
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development of the MTPV systems. 



III. CONCLUSION 
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FIG. 3: (Color online) q" , calculated for the vacuum gap 
(blue) and the gap filled with CNTs (red), (a) - d = 100 
nm, (b) - d = 1 /im. 



In this study we have shown the way to revolutionary 
increase the efficiency of microgap TPV systems inserting 
between the hot and the photovoltaic bodies a metama- 
terial which having the negligible thermal conductance 
converts the evanescent waves into propagating ones in a 
very broad range of spatial frequencies. We have shown 
that this conversion can hold over a wide frequency band 
and lead to a dramatic (2-3 orders) enhancement of the 
total thermal transfer across the micron gaps. 
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